XY model in small-world networks 
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The phase transition in the XY model on one-dimensional small- world networks is investigated by 
means of Monte-Carlo simulations. It is found that long-range order is present at finite temperatures, 
even for very small values of the rewiring probability, suggesting a finite-temperature transition for 
any nonzero rewiring probability. Nature of the phase transition is discussed in comparison with 
the globally-coupled XY model. 
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Examples of complex networks are abundant in many 
disciplines of science and have recently attracted much 
attention jjj. Among interesting phenomena is the so 
called small-world effect that a randomly chosen pair of 
nodes can be connected by a remarkably small number 
of intervening nodes. This effect, first noticed by Watts 
and Strogatz (WS) §, can be observed in a variety of 
real network systems such as the World-Wide Web, so- 
cial networks, and scientific collaboration networks P,H. 
The WS model is based on a locally highly connected 
regular network, in which some of the links are randomly 
"rewired" , creating long-range "shortcuts" . In such a 
model network the small- world effect is usually measured 
by the scaling behavior of the characteristic path length 
£, defined to be the average of the shortest distance be- 
tween two nodes: I ~ logiV with the network size (i.e., 
the number of nodes) N. Noteworthily, the small-world 
phenomena in the WS model emerge even at very small 
values of the rewiring probability P ~ 0(N~ 1 ), which 
implies that the global feature of the network is altered 
dramatically in the presence of only a tiny fraction of 
shortcuts. 

From the point of view of statistical physics, the above 
small- world effect can imply the emergence of global co- 
herence in the presence of shortcuts: In the absence of 
shortcuts, global coherence is difficult to achieve since 
the information to make each element have the same 
state should travel long distance of the order of the net- 
work size. As the number of shortcuts is increased, on 
the other hand, long-range connections become available, 
assisting the system to behave as a whole. The sig- 
nificance of such a shorter global length scale has been 
tested with many statistical physical problems studied on 
a small- world network topology. Included are the signal 
propagation speed j| , synchronizability || , dynamics of 
Hodgkin-Huxley neurons j^j , epidemiology |],|7j , percola- 
tion and — most relevant to the present study — the 
Ising model, where a nonvanishing order parameter has 
been demonstrated in the presence of a vanishingly small 
fraction of shortcuts jl(J . 

The XY model, describing two-dimensional (2D) spin 
(i.e., planar rotor) systems or superconductors and super- 
fluids, is one of the most well-known systems in statistical 



physics. In one and two spatial dimensions, the Mermin- 
Wagner theorem predicts that the XY model with only 
local interactions should not possess long-range order at 
any nonzero temperatures. Nevertheless, in two dimen- 
sions, it has been shown that although true long-range 
order is not possible, quasi-long-range order appears at fi- 
nite temperatures flli| . As the XY model includes more 
and more long-range interactions, one expects that the 
system can display true long-range order. In this work 
we study the XY model on the WS model of small-world 
networks and investigate the possibility of long-range or- 
der as the number of shortcuts or long-range interactions 
is increased. 

The Hamiltonian for the XY model reads 



cos 



;)> 



(1) 



where <pi is the angular direction of the 2D spin or the 
phase of the superconducting order parameter at node i. 
The coupling matrix J t j is given by 



J, if i and j are connected, 
0, otherwise. 



(2) 



For example, in the standard 2D XY model on a square 
lattice, we have = except for the four nearest neigh- 
bors of site i. In the WS model for the small- world net- 
work the regular network with N sites is first con- 
structed by connecting each site to its 2k nearest neigh- 
bors. Each local link is visited once, and then with the 
rewiring probability P it is removed and reconnected to a 
randomly chosen site. In this manner, a small-world net- 
work with size N is constructed for given k and P. We 
then study the phase transition of the XY model on this 
small- world network by measuring the order parameter 



(m) = 



N <V 



(3) 



where (• ■ •) denotes the thermal average and the aver- 
age over different network realizations, denoted by [■••], 
is also to be taken. When P = 0, the system behaves 



1 



as the one-dimensional (ID) XY model, displaying no 
long-range order at finite temperatures. For P ^ 0, the 
system has many long-range interactions in the thermo- 
dynamic limit and the Mermin- Wagner theorem cannot 
disprove the existence of long-range order. In particular, 
for P = 1 , it is on a random network with the number of 
randomly connected links being kN. 

We perform extensive Monte-Carlo (MC) simulations 
with the standard Metropolis local update algorithm, to 
measure various thermodynamic quantities and moments 
of the order parameter. For a thorough analysis of the 
MC results, we detect the phase transition in three differ- 
ent methods: First, Binder's fourth-order cumulant fill 



U N {T) = 1 



[(to 4 )] 
3[(to 2 )] 2 



(4) 



for different sizes N should have a unique crossing point 
at the transition temperature T c . Expanding the cumu- 
lant near T c , we write 



U N {T) «17* + E7i 1 



T 



AU N = Ujfin) - U N (T 2 ) tx N 1 /*, 



(5) 



(6) 



with T\ and T 2 (> Ti) picked near T c . From the above 
expression, the critical exponent v, describing the diver- 
gence of the correlation volume £y at T c [ [13||l4| ] 



IT - T r 



(7) 



can be determined. Our second method is based on the 
finite-size scaling of the order parameter, which exhibits 
the critical behavior [(m)] ~ (T c — T) 13 in the thermo- 
dynamic limit. In a finite-sized system, we expect the 
behavior to [(to)] = (T c - T)Pf{T,N) with a function / 
of at most two arguments T and N. We then use the 
standard finite-size scaling idea that the ratio of the cor- 
relation volume to the system size gives the argument 
of the scaling function, and obtain the finite-size scaling 
form 



[(to)] = N-P/°g ((T-TJA 1 / 5 



(8) 



which leads to a unique crossing point at T c in the plot 
of \{m)}N^/ 9 versus T. As discussed later in the present 
work, the phase transition on the small-world network 
is of the mean-field nature, which provides us the third 
method of locating the transition temperature: With the 
mean- field value of the critical exponent a = 0, the spe- 
cific heat neither diverges nor vanishes, remaining finite 
possibly with jump discontinuity. Accordingly, it is plau- 
sible to write the finite-size scaling form of the specific 
heat as 



C v = h((T-T c )N 1 ' 9 }, 



(9) 



which again crosses at T c . In the same way as in Eq. (0), 
the expansion of C v near T c can be used to determine v. 

All the above finite-size scaling forms are based on the 
assumption that we have only two spatial scales in the 
system: the size N and the correlation volume £y diverg- 
ing at T c . In fact it is known that the small- world net- 
work has an additional length scale, the typical distance 
between the ends of shortcuts, given by £ = (kP)^ 1 [0. 
Accordingly, in the presence of the three competing scales 
(iV, £, £), the standard finite-size scaling function should 
take the form x(£/Af, £/AT). Here we focus on suffi- 
ciently large systems with N much larger than £, where 
X(€/N>C/N) ma y be approximated as %(£/iV,0). This 
leads to the above-mentioned finite-size scaling forms 
without C [Eqs. (||), and @], which become more 
precise as the system size grows. This is manifested be- 
low when the MC results are presented. 

The small-world network, on which the XY model 
is defined, is characterized by the rewiring probability 
P and the local interaction range k. We focus only 
on the case k = 3; however, we believe that qualita- 
tively the same conclusion should hold for other values 
of k (> 1) p5[ . From the MC simulations, we have con- 
firmed that the network with P = does not exhibit 
long-range order at finite temperatures, in accordance 
with the Mermin- Wagner theorem. In the opposite case 
with P = 1, where all kN (= 37V) connections are long- 
ranged shortcuts, the system is found to undergo a well- 
defined finite-temperature transition. We first present 
our MC results for P = 0.2 in Fig. |l|. The finite-size scal- 
ing of Binder's cumulant, the order parameter, and the 
specific heat all reveals unanimously T c ks 2.235 in units 
of J/ks, confirming the presence of a finite-temperature 
phase transition in the XY model on the small-world 
network. In particular, the obtained critical exponents 
(3 = 1/2, V = 2, and a = establish the mean-field na- 
ture of the transition. Also confirmed is the validity of 
our finite-size scaling forms, with the additional length 
scale C reflected in the drift of the crossing points at small 
sizes. As the size N is increased, the finite-size effects as- 
sociated with £ reduce and the finite-size scaling forms 
based on the assumption of only two scales N and £y 
prevail, leading to well-defined crossing points at larger 
sizes. 

For comparison, we also study the globally coupled 
XY model, where Jij — J = 0(1/ N) for all i and j and 
T is measured in units of NJ/ks- We use the saddle- 
point method accompanied by the Hubbard-Stratonovich 
transformation fl6f| , and obtain various thermodynamic 
quantities such as the specific heat and the order pa- 
rameter in Eq. (||). We also perform the MC simulation 
and compare the results in Fig. |^, which exhibits excel- 
lent correspondence between the two. As expected, the 
globally-coupled XY model indeed possesses the mean- 
field transition, which is the same as that in the XY 
model on the small-world network. It should be noted 
that the number of connected links (i.e., the number 
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of nonzero elements in the interaction matrix J^) in 
the globally-coupled XY model is given by N(N—l)/2, 
which is in sharp contrast with the number O(N) in 
the small- world networks. This indicates that the XY 
model can be made to have global phase coherence (i.e., 
true long-range order) at finite temperatures by rewiring 
O(N) connections, with the total number of interactions 
still kept O(N) instead of <D(N 2 ) as in the globally- 
coupled XY model. 

What is the minimum value of P to have such long- 
rage order at finite temperatures? To answer this ques- 
tion, we have also performed MC simulations at various 
values of P. As the rewiring probability P is reduced, 
the finite-size effects due to the length scale C grow and 
become no more negligible, making it necessary to per- 
form simulations on larger systems. Since it is observed 
that the finite-size effects are much more discernible in 
Binder's cumulant and in the specific heat than in the 
order parameter, we estimate T c only from the crossing 
point in the plot of [(m)]^ 1 / 4 versus T. For P > 0.03, 
we have confirmed that as the system size grows, all the 
three methods to determine T c give the same value. As 
an example, we show in Fig. |^ the determination of T c 
for P = 0.05. 

Finally, Fig. | summarizes the results and presents the 
phase diagram for the XY model on small-world net- 
works with k — 3. We observe that the phase boundary 
is very well described by the logarithmic form: T C (P) « 
0.41 InP + 2.89, as shown in the inset of Fig. f|, without 
an obvious reason. A naive extrapolation of this form 
predicts a very small value of the critical rewiring prob- 
ability: P c w 0.001. It is, however, likely that there 
takes place deviation from the logarithmic dependence 
at small P and any nonzero rewiring probability (P > 0) 
presumably supports long-range order at sufficiently low 
but finite temperatures, resulting in P c = 0. One may 
provide a simple argument in favor of this expectation: 
Since the spins separated within the correlation length 
are correlated with each other, the presence of shortcuts 
in such a way that the typical distance £ between the ends 
of shortcuts is smaller than the correlation length £ (of 
the ID system) does not affect the system substantially, 
leaving the ID nature intact. When £ grows beyond the 
correlation length, on the other hand, the long-range in- 
teractions via shortcuts come into play, giving rise to the 
mean-field character. Accordingly, the crossover between 
the ID behavior and the mean-field one is given by the 
condition £ ss £, which also describes the so-called small- 
world transition between the two geometrical regimes, 
the large- world regime and the small- world one JTtJ . Rec- 
ognizing the low-temperature behavior of the correlation 
length £ ~ T _1 jl8| in one dimension [P — 0) and the 
typical distance £ ~ P~ x , we thus expect the behavior of 
the transition temperature T c ~ P in the limit of small 
P |ll| . Here it is of interest to note the difference from 
the Ising model on small-world networks. The latter ex- 
hibits the behavior T c ~ — (InP) -1 [jl0|, which originates 
from the exponential temperature dependence in the ID 



Ising model, £ ~ e a / T with the constant a proportional 
to the coupling strength. 

In conclusion, we have studied the XY model on small- 
world networks through the use of the standard Monte- 
Carlo simulation method. At all nonzero values of the 
rewiring probability P considered in this work, the finite- 
temperature transitions have been detected by Binder's 
cumulant, the specific heat, and the order parameter, 
together with appropriate finite-size scaling forms. The 
phase transition has been shown to be described by the 
mean-field critical exponents, (3 = 1/2, D = 2, and a = 0, 
irrespective of the value of P. The existence of the addi- 
tional length scale corresponding to the typical distance 
between the ends of shortcuts is manifested by the drift 
of crossing points at smaller sizes. The question as to 
the value of P c , below which long-range order does not 
emerge at any finite temperature, still remains to be an- 
swered. The phase diagram in Fig. ^ provides a very 
small upper bound, P c < O(10 -3 ), which is likely to sug- 
gest P c = 0. In this regard, it is revealing that the small- 
world transition, where the characteristic path length 
£ undergoes a change of behavior, occurs at P = in 
the thermodynamic limit |l7lj . We thus suspect that the 
small-world transition and the order-disorder transition 
are intimately related, detailed investigation of which is 
left for further study. 
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FIG. 1. The XY model on the small- world network with 
P = 0.2 and k = 3 displays a mean-field phase transition at 
the finite temperature T c : (a) Binder's cumulant Un has a 
unique crossing point at T c w 2.235 (in units of J/ks) and 
the critical exponent v ~ 2.0 is obtained from Eq. (^) (see 
the inset), (b) Finite-size scaling of the order parameter [see 
Eq. ||] again yields (3 « 0.5 and v » 2.0 with T c « 2.235. 
(c) Specific heat C v (in units of ks) also has a crossing point 
at T c « 2.235, as implied by Eq. (|). 
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FIG. 2. Specific heat C v (in units of feg) and the order pa- 
rameter (m) of the globally-coupled XY model. Compared 
are C v and (m) from analytic calculation (full and dotted 
Unes, respectively) and from Monte-Carlo simulations for the 
system size N = 6400 (empty squares and filled circles). The 
phase transition is of the mean-field type characterized by 
(rn) ~ (T c - T) p with ,3 = 1/2 and the jump in C v at the 
transition (T c = 1/2 in units of JN/ks)- 
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FIG. 3. XY model on the small- world network for 
P = 0.05. From the crossing points in the scaling forms of 
(a) the order parameter \(m)\ and (b) the specific heat C v (in 
units of fcs), T c ~ 1.66 (in units of J/ks) is obtained. The 
drift of the crossing point at small sizes is larger in the specific 
heat (b) than in the order parameter (a). 
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FIG. 4. The phase diagram of the XY model on the 
small- world network with k — 3, drawn on the plane of the 
rewiring probability P and the temperature T (in units of 
J/fcs). The data points have been obtained from the fi- 
nite-size scaling form of the order parameter in Eq. (^) and 
the error bars are smaller than the size of the symbol. The 
solid line, which is only a guide to eyes, gives the boundary 
separating the disordered high-temperature phase (D) from 
the low-temperature phase with the non- vanishing order pa- 
rameter (O). Inset: The phase boundary is well described by 
the form T ^ 0.41 InP + 2.89 except for the point P = 1. 
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